Probing Anomalous Longitudinal Fluctuations of the Interacting Bose Gas via 

Bose-Einstein Condensation of Magnons 
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The emergence of a finite staggered magnetization in quantum Heisenberg antiferromagnets sub- 
ject to a uniform magnetic field can be viewed as Bose-Einstein condensation of magnons. Using 
non-perturbative results for the infrared behavior of the interacting Bose gas, we present exact re- 
sults for the staggered spin-spin correlation functions of quantum antiferromagnets in a magnetic 
field at zero temperature. In particular, we show that in dimensions I < D < 3 the longitudinal 
dynamic structure factor S«(q,u) describing staggered spin fluctuations in the direction of the stag- 
gered magnetization exhibits a critical continuum whose weight can be controlled experimentally by 
varying the magnetic field. 
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Because magnons in ordered Heisenberg magnets are 
bosonic quasi-particles, it is natural to expect that un- 
der certain conditions these systems can be used to study 
general properties of interacting Bose gases, such as the 
phenomenon of Bose-Einstein condensation (BEC). In 
fact, the renewed interest in BEC in recent years has 
also led to considerable activity in the field of magnon 
BEC which we abbreviate MBEC below. 

Experimentally, MBEC has been observed in two dif- 
ferent classes of systems. The first are spin-singlet sys- 
tems [l| such as TICUCI3 or Haldane gap spin chains, 
where MBEC corresponds to the collapse of the singlet- 
triplet gap at a critical magnetic field and the emergence 
of long-range magnetic order for larger fields. The second 
class of systems, which we shall consider in this work, 
are easy-plane quantum antiferromagnets (QAFs) such 
as CS2CUCI4 subject to a magnetic field h perpendicular 
to the easy plane 0. If h exceeds some critical field h c , 
the ground state is a saturated ferromagnet. For h < h c 
the [/(l)-symmetry of the Hamiltonian is spontaneously 
broken and antiferromagnetic long-range order emerges. 
At zero temperature, the boson density vanishes at the 
phase transition so that many-body interactions are not 
relevant at the zero temperature quantum critical point. 
However, the finite temperature transition belongs to the 
same universality class as BEC in the interacting Bose 
gas. One advantage of studying BEC in magnetic sys- 
tems is that the magnetic analogon h c — h of the chemi- 
cal potential can easily be controlled experimentally. An- 
other advantage, which has attracted little attention to 
date, is that in MBEC the phase of the condensate has 
a direct physical interpretation as the orientation of the 
staggered magnetization within the easy plane. 

The concept of MBEC has been formulated theoret- 
ically many years ago [1, 0, H, @|, where the focus was 
mainly on the nature of the quantum phase transition. 
However, even away from the critical point, perturba- 
tion theory for the interacting Bose gas is plagued by di- 



vergences originating from the gapless dispersion of the 
Goldstone mode, giving rise to non-analytic behavior of 
some boson correlation functions 0, H, HJ. Yet, in all 
physical observables of the Bose gas these divergences 
cancel, so that they seem to be little more than a math- 
ematical curiosity. Here, we point out that, in contrast 
to the Bose gas, in MBEC the anomalous behavior aris- 
ing from the divergences can be directly observed exper- 
imentally. We shall exploit recent renormalization group 
results 7] to obtain the infrared behavior of the stag- 
gered spin-spin correlation functions of QAFs in a uni- 
form magnetic field at zero temperature. In particular, 
we shall show that in dimensions 1 < D < 3 the lon- 
gitudinal part Sit (q, lu) of the dynamic structure factor 
(corresponding to spin-fluctuations parallel to the stag- 
gered magnetization) exhibits a critical continuum for 
small wave-vectors \q\ and frequencies ui, which can be 
measured via neutron scattering. 

We start from the Hamiltonian of the spin-S" QAF on a 
D-dimensional hypercubic lattice with N sites and lattice 
spacing a in a uniform magnetic field in z-direction, 



(1) 



where the sums are over all sites of the lattice, and Si 
are spin operators with Sf — S(S + 1). We assume 
nearest neighbor coupling, i.e., — J > if i and j 
are nearest neighbors, and Jy = otherwise. For suf- 
ficiently large h, the ground state of Eq. Q is a sat- 
urated ferromagnet with magnetization parallel to the 
field. To obtain the low-lying magnon excitations, we ex- 
press the spin operators in terms of canonical boson op- 
erators bi using the Holstein-Primakoff transformation, 
Sr = V2Sb\[l - n i /{2S)} 1 ' 2 = (S+)t, S? = S - n h 
where u, = fej&j. Expanding the square roots and ne- 
glecting terms involving six and more boson operators 
we obtain H w Eq + H2 + H4. The constant part is 
E a = N[J S 2 /2 - hS], where J = J fe=0 with J k = 



2 



J2j e j Jij . The quadratic and the quartic contribu- 
tions to the Hamiltonian are 

H 2 = --^2 Jij [ni + rij - b\bj - b]bi] + h^m,{2) 



Hi = ~ ^2 Jij&niUj ~ b\riibj - b^n^] 



(3) 



Anticipating antiferromagnetic symmetry breaking, it 
is convenient to define the Fourier transform in the 
reduced (antiferromagnetic) Brillouin zone with two 
branches of magnon operators bka labelled by momenta 
fc in the reduced zone and a = ±, 



b ka = N-v*[y2e 



ieB 



(4) 



Here i £ A/B denotes a summation over the sites of only 
one of the sublattices A or B. The quadratic magnon 
Hamiltonian is then diagonal, H 2 = J2ka( e k<? ~ vWhtfiko, 
with magnon dispersion tka = (Jo + uJk)S and /i = 
h c — h, where h c — 2JqS. For small wave- vectors tk,- ~ 
fc 2 /(2m), with effective mass m = (2JSa 2 )^ 1 , while the 
symmetric mode ek,+ is gapped, tk=a,+ = h Cl so that 
£fc=o,+ — fJ> = h. Note that the gapless mode a = — 1 
describes staggered spin fluctuations. 

Since in this work we are interested in the infrared be- 
havior of correlation functions at energy scales smaller 
than h, we shall ignore the gapped mode. To derive 
the interaction part of the Hamiltonian of the gapless 
magnons, we substitute Eq. ((4| into Eq. ((3]) and neglect 
all terms involving the operators bk,+ - Defining the field 
operators tpk — V 1 / 2 bk,-, where V — Na D is the volume 
of the system, we obtain for N — > oo, 
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(5) 



where f. = j d D k/(2n) D and the effective interaction is 
U q = 6(A — |q|)xo 1 - Here A w a -1 is an ultraviolet 
cutoff and xo = (2Jo a£> ) _1 is the leading large-5 1 result 
for the uniform transverse susceptibility in x-direction, 
which is perpendicular to the staggered magnetization 
and to the magnetic field. This completes our mapping 
of the spin system onto an interacting Bose gas. 

For /j > the field operator ipk=o acquires a vac- 
uum expectation value, corresponding to the emergence 
of long-range antiferromagnetic order. The spin config- 
uration in the ground state is then canted, as shown in 
Fig. HJ One should keep in mind, however, that Eq. (JSJ) is 
accurate only for small 6, since otherwise the Bose gas is 
no longer dilute such that a three-body interaction must 
be added. However, the qualitative behavior of corre- 
lation functions is largely determined by symmetry Q, 
so that our results derived below are expected to remain 




A B ---^j A B^\, 

FIG. 1: (Color online) Spin configuration in the ground state 
of a QAF in a uniform magnetic field for h < h c . The local 
moments (Si) are represented by dashed arrows. We choose 
the coordinate system such that the magnetic field h points 
in the ^-direction and the staggered magnetization M s points 
in the y-direction. The sublattices are labelled A and B. 



valid also for larger canting angles. Of course, if 9 is 
not small it is better to quantize the spin operators in 
a suitably defined tilted basis U|- For simplicity, 



we shall focus here on the regime 9 <^ 1 where our ap- 
proach based on the Holstein-Primakoff transformation 
with quantization axis along the direction of the magnetic 
field is thus quantitatively accurate for large S. Alterna- 
tively, we could represent the spin operators in terms of 
projected bosons which account for all l/S cor- 

rections at small 9. However, in the spirit of Ref. [l3j ]. 
we take the finite corrections due to higher orders in 1/5 
implicitly into account by expressing our final result in 
terms of the true spin-wave velocity c and the spin sus- 
ceptibility x, which can be obtained from experiments. 

To leading order in l/S, the classical tilt angle 9 is eas- 
ily obtained by substituting = (27r)- D (5(fc)'0o + 
in Eq. (J5]) and demanding that the coefficient of the 
term linear in A-0fc vanishes. This yields the well known 
14| relation IV'opfo = To translate this into a 
condition for 9, we note from Fig. [1] that for large S 
the tilt angle is related to the staggered magnetization 
M s = V- 1 J2i Ci(Si) via 9 « M s /s. Here s = S/a D is the 
spin-density of the saturated magnetic state, and Q = 1 
for i G A and Q = -1 for i G B. Since M 2 = 2s|-0o| 2 to 
leading order in spin-wave theory, we obtain for large S 



Mils 2 » 2p /s = 2(1 - h/h c ) , 



(6) 



where po = \4>o\ 2 is the condensate density. Eq. ([6]) 
agrees with the small-f? expansion of the classical re- 
sult cos0 = h/h c [HI]. To obtain the magnon spectrum 
for large S, we neglect all terms higher than quadratic 
in the Aipk- After a Bogoliubov transformation, one 
finds that for h < h c the magnon spectrum is deter- 
mined by ejk = (|^r) + Cgfc 2 which for small fc yields 
(k = co\k\ + C(fc 3 ), with the magnetic field dependent 
spin-wave velocity cq = yjijm = 2\f~D3Sa9. 

We now focus on the spin-wave interactions, which are 
relevant in D < 3. The substitution iph = (27r)- D (5(fc)'0o + 
Aipk generates various terms cubic and quartic in 
which are rather tedious to work with. A perturbative 
treatment of these terms is plagued by infrared diver- 



3 



gences. However, Ward-identities associated with the 
{/(l)-symmetry of the Hamiltonian can be used to ob- 
tain the infrared behavior of some correlation functions 
without resorting to perturbation theory. As pointed out 
in Rcf. 7], the constraints imposed on the correlation 
functions by the Ward identities are more transparent if 
one expresses ?pk in terms of two real fields, correspond- 
ing to transverse and longitudinal fluctuations. Recently 
we have proposed a similar parameterization of the spin- 
wave expansion in Q AFs [HI, [l(| , which amounts to ex- 
pressing the field operator tpf. in terms of two canonicall 
conjugate hermitian operators 11^ and as follows 
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(7) 



Both fields represent staggered spin fluctuations in the di- 
rections perpendicular to the magnetic field, where life is 
transverse and longitudinal with respect to the stag- 
gered magnetization. More precisely, in the coordinate 
system shown in Fig. [T]we have for 5* — > oo, 



— Yb< 



— ik-'- 



(8) 

This equation becomes quantitatively accurate even for 
finite 5* in the dilute limit h — > h c which can be shown 
using projected boson operators 12J. The condensed 



phase corresponds to antifcrromagnetic order, ($fc) = 
(2t:) d 5 (k)<j) , where O = M 2 /s for large S. Note that 
the total density of the Bose gas corresponds in the un- 
derlying spin system to p = s— V^ 1 J^iv^f)) which van- 
ishes for h > h c and satisfies p <C s if h is only slightly 
smaller than h c . Corrections to Eq. |[8j would be asso- 
ciated with terms which are of linear or higher order in 
the boson density and become negligible for h — > h c . 

The infrared behavior of correlation functions is most 
conveniently derived within a functional integral ap- 
proach [3, EH ■ The system is then described in terms 
of an action 5 [II, $], which is a functional of c- number 
fields depending on imaginary time r. Formally, these 
fields can be obtained from the field operators by sub- 
stituting f[ fe -> n fc (r) and $ fc -»■ (2ir) D 6(k)(f> Q + $ fc (r). 
Introducing the Fourier transforms in frequency space, 
U K = J dre^ r n fc (r) and § K = J dre MT $fc(r), the 
Gaussian part of the action can be written as 



1 



K 



5 [n,$] = - / X oc 2 k 2 U_ K U K + Xo^-k^k 



(9) 



Here K = (fc, iu>) is a collective label for momenta and 
frequencies, and J K — f J k - The correlation functions 



are within Gaussian approximation given by 



(RkHk')o = Sk,-k' 



(ILk®k>)o = $k-k> 



Xo 1 



uS 1 + Cpfc 2 



uS 1 + Cpfc 2 



(10a) 
(10b) 

(10c) 



where 6 k ,-k> = (2tt) d+1 S(uj + u')S(k + k'). It turns 
out that the Gaussian approximation is qualitatively cor- 
rect for the transverse correlation function (TIk^Ik 1 ) and 
for the mixed correlation function (Hk&k 1 ), where (...) 
denotes the full thermal average. Using the results of 
Ref. 0] , we obtain for the true infrared behavior of these 
correlation functions 

-l 



(n x n. 



X 



(Uk^k') = 5k,-k> 



jj* + c 2 k 2 
Zulu 



■z 2 k 2 1 



(11a) 
(lib) 



where c is the renormalized spin-wave velocity, \ — 
M 2 p/(2smc 2 po) = x^ZpjZ 2 is the renormalized trans- 
verse susceptibility, and po and p are the true condensate 
density and total density. For convenience we have intro- 
duced the dimensionless factors Z p = p/po, Z c = c/cq, 
and Z\\ = (mc 2 / p)dpo/dp 1 all of which approach unity 
in the classical limit S — ► oo. We emphasize that in di- 
mensions D < 3 some Feynman diagrams contributing 
to the above correlation functions are infrared divergent. 
However, the Ward identities guarantee that all di- 
vergences cancel, so that the only difference between the 
exact results (lllai II lb[) and the corresponding approx- 



imate expressions (llOal I10b[) are finite renormalization 
factors. The exact results (jllal lllbp depend only on 
parameters which can be expressed in terms of thermo- 
dynamic derivatives. 

The important point is now that the Ward identities 
do not impose similar constraints on the longitudinal cor- 
relation function ($k&k>), which in D < 3 is dominated 
by a non-analytic term arising from the re-summation of 
the leading infrared divergences. Using the results de- 
rived in Ref. Q we obtain r 



-K D+ 




1 Z pP0 



(12) 



where K D = 2 1 - d tt- d / 2 /T(D /2). The non-analytic con- 
tribution to the longitudinal correlation function of the 
interacting Bose gas has first been discussed by Weich- 
man and later in Ref. The field 4>k has thus a 
finite anomalous dimension and its effective action can- 
not be approximated by a Gaussian. This reflects the 
general behavior of systems with broken continuous sym- 
metries where Goldstone modes lead to anomalous lon- 
gitudinal fluctuations 0, [H, S HH, HI . Note that the 



4 



magnetic field dependent microscopic momentum scale 
mc is for large S approximately mc ~ \J~D9ja. Using 
a generalized Ginzburg criterium Q, we find that the 
non-analytic corrections in Eq. (|12|) become important 
for |fe| < k G , where k G w mc [(mc) D / po] 3 ~ D for D < 3, 
and kc w mcexp[— po/(wc) 3 ] for £) = 3. 

Even though we have derived Eqs. (|lla[|llb| and (fT2")) 
only for small canting angle 0, they remain valid even for 
moderate values of 9, as long as three-body interactions 
are negligible. Using Eqs. (|8|12p . we find for the longi- 
tudinal staggered structure factor for 1 < D < 3 and 
< uj/c < k G , 



Cr 



7 2 

£lL 

2 

(mc) 3 



c|fc|*(w - c\k\) 
G(o;-c|fe|) 



(13) 



where C D = ifr,+i[7r(3 - D)/2}- 1 sin[7r(3 - L>)/2]. In 
particular, C 3 = 4 = (87T 2 )- 1 and C 2 = tt" 3 . Eq. 112) 
is the main result of this work. For small |fe| < k G 
the critical continuum represented by the last term in 
Eq. (fT3"|) carries most of the spectral weight: Denot- 
ing by I c the contribution from the continuum part in 
Eq. (fT3"|) to the energy integrated spectral weight, and by 
Is the corresponding contribution due to the <5-function, 
we find I c /I s oc (k G /\k\) ln(fc G /|fc|) for D = 2 and 
oc (fcc?/|fc|)(mc) 3 //9o for D = 3. In both cases 
the continuum part dominates for |fe| < kc- However, 
in three dimensions k G oc 9 exp[— S/(6>/36)] is exponen- 
tially small for h w h c so that the region |fe| < k G is 
difficult to resolve experimentally. On the other hand, 
in D = 2 we have k G w 4V2<9/S a so that in this case 
the critical continuum should be accessible via neutron 
scattering at scales |fc| < k G . 

In summary, we have shown that in dimensions 1 < 
D < 3 the longitudinal dynamic structure factor of 
QAFs subject to a uniform magnetic field exhibits a field- 
dependent critical continuum as described by Eq. (|13p 
which is closely related to the anomalously large longitu- 
dinal fluctuations in the condensed phase of the interact- 
ing Bose gas. While in the latter case these fluctuations 
cannot be directly measured, the magnetic realization of 
BEC offers a direct experimental access to these fluctu- 
ations. Both in D = 2 and D = 3 the critical contin- 
uum dominates the energy integrated longitudinal struc- 
ture factor if |fe| is smaller than the Ginzburg scale k G . 
At least in two-dimensional QAFs the critical continuum 
should be observable in neutron scattering experiments. 
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